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Abstract. The measurement of high-dimensional entangled states of orbital angular momentum prepared 
by spontaneous parametric down-conversion can be considered in two separate stages: a generation stage 
and a detection stage. Given a certain number of generated modes, the number of measured modes is 
determined by the measurement apparatus. We derive a simple relationship between the generation and 
detection parameters and the number of measured entangled modes. 
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1 Introduction 

Entangled states are a distinctive feature of quantum me- 
chanics. Their use can lead to important technological ad- 
vances in communication, security and, ultimately, com- 
puting [T]. Entanglement in a high-dimensional Hilbert 
space means a high effective number of entangled modes 
that can be used to achieve a high shared information 2 . 
It is therefore of great importance to choose the proper 
basis in which to detect entangled modes. The Schmidt 
basis is that which yields the maximum shared informa- 
tion [3]. For most entangled states, however, the elements 
of the Schmidt basis cannot readily be measured, perhaps 
because the size of the components in the detection appa- 
ratus do not match any possible detection mode 0. 

It is well-known that light can carry orbital angular 
momentum (0AM) and that this property is associated 
with a helical phase front ^ (and papers reprinted therein). 
Optical modes carrying 0AM include the Laguerre-Gaussian 
modes [7] and also the Bessel beams. Of central interest 
to us in this paper is the fact that photon pairs produced 
in spontaneous parametric down-conversion (SPDC) are 
naturally entangled in their 0AM |8I9I10) . One clear man- 
ifestation of this entanglement is the existence of an EPR 
paradox for the 0AM and its conjugate quantum variable, 
the azimuthal angle [12113115] . 

The 0AM is conserved in the down-conversion process 
and hence for a Gaussian {£ = 0) pump, the 0AM of the 
signal and idler fields are perfectly anticorrelated. There 
are also correlations on the radial direction (as quanti- 
fied, for the Laguerre-Gaussian modes, by a radial index 
p) [IS] but these will not concern us in this paper. Our 
central concern will be the number of entangled lowest 
order (p = 0) Laguerre-Gaussian modes generated in a 
down-conversion experiment. The typical setup that we 
consider is a type-I or type-II, degenerate SPDC setup. 



We work in the regime of undepleted pump and we ne- 
glect eventual anisotropics of the down-converted beams. 

We find that, for any given set of generation parame- 
ters (pump waist Wp, wavelength A, crystal length L) the 
detection apparatus can be prepared in a way that max- 
imises the measured number of entangled modes and that 
two important parameters are 7, the ratio of the width of 
the pump beam to the width of the detection modes, and 
Lfl, the length of the crystal normalised to the Rayleigh 
range of the pump beam: 



AT ^ 

— — and Ln = — . 

Ws,i Zr 



(1) 



Where the Rayleigh range is zr = ~x^- In this paper we 
assume that the signal and idler modes have the same 
width so that Ws = Wi and 7s = 7i = 7. 

The precise calculation of Ws.i depends upon the de- 
tails of the detection system. Our analysis can be applied 
if the back-projected detection mode size, Wi,s, is approx- 
imately ^- independent over the range of 0AM of interest, 
and if the modes with p ^ couple only weakly with the 
fundamental mode of the fibre that carries the signal to 
the coincidence counter. 

We investigate the Lr dependence of the 0AM band- 
width, while recognising that many experiments operate in 
a regime where Lr <^ 1 |8ll7l ?i^.?]. In the short crystal 
limit and near to collinearity the familiar sine phase can be 
dropped [H]. One can then obtain an analytical form for 
the down-converted state |16I22] and its extension to non- 
Gaussian pump beams |23| . Our aim in this paper is to 
go beyond these existing analyses and to explore regimes 
in which the sine phase matching term becomes signifi- 
cant, which leads to the exact analytical expression ^ 
and to the characterisation of the detection parameters. 
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We present both an analytical treatment and also a simple 
geometrical argument for our results. 

The second section of the paper specifies the definitions 
of the various bandwidths which are used. The third sec- 
tion contains the analytical approach to calculate the pro- 
jection amplitudes. The fourth section contains the geo- 
metrical approach to calculate a simple formula that gives 
the measurement bandwidth. The fifth section contains 
the interpretation of the results and the conclusions. 



2 Definition of bandwidths 

For a distribution of probabilities, in our case for the OAM 
of the signal or idler photon in SPDC, we can define a 
number of statistical measures. For high-dimensional en- 
tanglement we require as many modes as possible to con- 
tribute to the state and, moreover, for these to contribute 
strongly, that is to have a significant probability. A simple 
and convenient measure of this quantity is the Schmidt 
number [M05] : 

K{{Pr}) ■■= ^ , (2) 

where the probabilities {pi} are, in our case, those for 
each of the OAM modes. The measure K gives the effec- 
tive number of contributing modes and hence the effective 
dimensionality of the system. In experiments, it is typical 
to quote the full- width at half maximum as the measure of 
the bandwidth (FWHM) so as to include only modes that 
are well above the noise fioor. FWHM should not be con- 
fused with K. For simple, symmetrical and single-peaked 
probability distributions, the Schmidt number provides a 
convenient measure of the bandwidth. The precise rela- 
tionship between the Schmidt number and the FWHM 
depends upon the detailed shape of the distribution but 
typical of our systems is that the K exceeds the FWHM, 
see figure [TJ For a distribution like this we can define an 
effective range of modes contributing to the state ranging 

from £i„ax to £min = "^max SUCh that K =1 + 2|Cax|- 

The generation bandwidth is the effective number of 
entangled modes generated in the SPDC process. As it 
does not depend on the detection apparatus, it is a func- 
tion only of the crystal length and of the size of the pump 
beam, combined into the quantity Ln, defined in eq. ([1]). 
This bandwidth can be thought of as the dimensionality of 
the entanglement in OAM and can be calculated through 
the Schmidt decomposition of the SPDC state [5]. More 
on the generation bandwidth is detailed in its derivation, 
in section 4.1. 

The measurement bandwidth represents the number of 
modes that a detector will measure in an experiment and 
depends on both the generated modes and on the overlap 
of these with the detection modes. In doing so, we need 
to consider the optics used to image the light onto the 
detectors and any restriction arising from this, such as a 
restriction to p ~ Laguerre-Gaussian modes. The over- 
lap between the generated modes and the back-projected 



detection modes needs to be maintained both in the im- 
age plane and in the far field plane of the crystal: a setup 
with high overlap in the image plane may still suffer from 
low overlap in the far field or vice versa and this would 
translate into a decreased modal sensitivity. This overlap 
requirement has a central role in the derivation of eq. (jlOp , 
which is based on the argument that the angular spread 
of a generated mode cannot exceed the natural spread of 
the down-conversion cone. In the next sections we will de- 
fine an image plane bandwidth and a far field bandwidth 
and, as we shall show, there is a natural way of combining 
the two. This geometrical result is strongly supported by 
the more complicated analytic result, which we evaluate 
numerically for a comparison in figure 3. 

3 Analytical treatment 

A direct calculation of the measurement bandwidth needs 
to consider the overlap between the SPDC state and a 
pair of joint detection modes |10I16| . This yields a series 
of complex measurement amplitudes {Ci} where £ labels 
each value of the OAM that was measured. The mea- 
sured Schmidt number (or the measurement bandwidth) 
is therefore given by the measure K applied to the set of 
projection probabilities 

K{{Pi}), where = (3) 

We seek to evaluate this quantity for a Gaussian pump 
laser, taking full account of the sine phase-matching term. 
In this way we extend the regime of validity of earlier 
calculations. 

We consider the measurement modes for the signal 
and idler fields to be a pair of Laguerre-Gaussian modes. 
The LG modes are characterised by two integers £ and p 
and a real positive number w, which represent the OAM 
quantum number, the radial quantum number and the 
Gaussian modal width, respectively. For simplicity, we set 
p = 0, which limits our analysis to modes with a single 
bright ring in the transverse plane. Many of our experi- 
ments are designed to detect p — modes with a higher 
efficiency, moreover, than higher-order modes. We note 
however, that modes with non-zero p are produced in the 
SPDC process [TB] and, indeed, it is these that makes it 
possible to observe entanglement of three-dimensional vor- 
tex knots in SPDC [25]. 

The SPDC wave function ip{cii, Qs), in momentum space, 
is written in the following way, where the subscripts s and 
i refer to signal and idler modes |10j : 

V'(q„q.)=iVe-^lq.+q=l%inc(^A|q,_q^|2^ . (4) 

Here q is the transverse component of the momentum vec- 
tor k, Wp is the pump width, L is the crystal thickness, kp 
is the wave vector of the pump. The first term corresponds 
to the transverse wavevector components of the pump, 
while the second term represents the phase-matching im- 
posed on the down-conversion process by the nonlinear 
crystal. 
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We consider each detection mode to be an LG mode 
with radial quantum number p — 0. In polar coordinates 
(p, if) in momentum space it has the form 



LGi{p,(p) = 



27r . 



pw 

71 



(5) 



The projection amplitude is therefore calculated by eval- 
uating the overlap integral of tp with two LG modes of 
opposite 0AM (because of angular momentum conserva- 
tion) [8l9ll0j . The result is found to be 



Lr 



■272 



Kl 



(6) 



We note that the first term in brackets corresponds to 
that obtained previously [l6i, specialised to equal signal 



and idler widths and p — modes. Here the function <l> 



Lr,1 



is the Lerch transcendent function of order (1, \£\ 
argument —27^^ }27j : 



^-(-272^, 1,1^1 + 1), c 



1) and 



(7) 



Note that ^^lior Lr = 0. 

Once Lji and 7 are specified, the amplitudes C^"-'^ 
are to be used in eq. (O, in order to calculate the mea- 
surement bandwidth. The dependence of the projection 
amplitudes on a transcendent function makes further an- 
alytical calculation difficult, and a numerical approach has 
to be employed. However, as the tails of the distribution 
of projection probabilities have a slow decay and there- 
fore an effect on the width even at high |^|, the numerical 
approach is slow, if an accurate result is sought. 



FWHM 




Fig. 1. An example of a distribution of \C^"''^\'^ for 
Lr = 0.001 and 7 = 2, obtained by calculating numer- 
ically the projection amplitudes between ^ = — 20 and 
£ = 20. The FWHM and the measurement bandwidth K 
are shown in blue and red, respectively. Note that K ex- 
ceeds the FWHM by about 2.5 times, giving an effective 
mode number of about 17 in this case. 



In figure [T] we give the probabilities for the angular 
momentum values ^ for Lr — 0.001 and 7 = 2. In this 



parameter range existing analytical expressions provide 
an excellent approximation |10I16| . 



4 Geometrical argument 

In this section we find an upper (and therefore lower) 
bound for the generated 0AM values, and for the mea- 
sured 0AM values. The measurement bandwidth that we 
calculate from such bounds matches the analytic result of 
the previous section and therefore allows to avoid calcu- 
lating numerically the distribution of projection probabil- 
ities. 




Fig. 2. The relation between a and Akz sets a natural 
upper bound to a for near-collinear emission. 



The phase-matching efficiency of the down-conversion 
process depends upon the axial mismatch Akz in wave 
vectors of the pump, signal and idler fields, and it is given 
by sinc^ {LAkz/2). When optimised for degenerate, near- 
collinear phase-matching, the signal and idler output is ob- 
tained over a narrow range of angles, a, for which LAkz < 
IT. With reference to figure [U for small a (which corre- 
sponds to being near to collinearity) we can write 



AK, 



a ko 



(8) 



It follows, therefore, that the allowed values of a are bounded 
from above: 



v2 < 



27r 



(9) 



For Laguerre-Gaussian modes, in the paraxial regime, 
we can define an effective local wavevector associated with 
the gradient of the phase. The helical form of the wave- 
fronts gives rise to an angular spreading of these such that 
at a distance r from the mode axis, the angular spread 
is P ~ t/kr [28 , which can be interpreted as the lo- 
cal spreading angle from the optical axis. The natural 
restriction on a imposed by the phase matching there- 
fore sets a limit /3 < a on the efficiency of production of 
the 0AM carrying beams, imposing a restriction on the 
generated 0AM bandwidth. Such restriction is a natural 
consequence of the fact that a generated mode cannot be 
more divergent than the down-conversion cone. The rela- 
tion /3 ^ a, using the definitions and bounds given above 
for P and a, can be rewritten as 



< 



Yl 



(10) 
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where we have made the approximation that fc^.i ~ kp/2. 
This relation is the starting point to calculate the genera- 
tion bandwidth and for the analysis in the far field of the 
image plane of the crystal. 



4.1 Generation bandwidth 

The beam size can be no bigger than that of the pump 
beam, i.e. r < Wp. Applying this bound to eq. (llOp we 
obtain an upper bound for the generated OAM value: 



(11) 



It follows, therefore, that the generation bandwidth is 




(12) 



This number represents the effective number of entangled 
OAM modes generated by the source obtained by remov- 
ing the p — restriction (as we are applying such restric- 
tion only to the measurement bandwidth) . Equivalently, it 
can be thought of as the bandwidth obtained by removing 
the restriction on 7, i.e. if one does detect p = 0, but with 
any 7. This way of thinking about Kg^n can be helpful, as 
it relates to a measurement scheme. The relation between 
-ftTgen and the total Schmidt number K or its azimuthal 
part ifaz [?] is not straightforward, because Kgen can be 
thought of in terms of a measurement with any value of 
7- 



4.2 Image plane bandwidth 

As anticipated in section 2, to calculate the measurement 
bandwidth we need to consider the overlap of the gener- 
ated field with the detection modes in the image plane 
of the crystal and in its far field. Intuitively, a detection 
system which has a good overlap in the image plane, but 
that detects light that only comes from a narrow spread of 
directions would restrict the measured bandwidth. A simi- 
lar restriction would occur for one that has a good overlap 
with the typical incoming angles of LG beams, but that 
has a poor overlap with the intensity in the image plane. 
It is clear that in order to optimise a detection system, 
both these quantities have to be taken into account. 

To calculate the overlap in the image plane it suffices 
to note that a, p = Laguerre-Gaussian mode with OAM 
number i and width w has its maximum intensity at a 
radius 



(13) 



For efficient conversion of pump to signal and idler we 
require that the pump, single and idler beams should all 
overlap, giving a restriction on the maximum size of the 
down-converted beams (r^^i < Wp) and hence an upper 



bound to the value of OAM in the plane of the crystal 
corresponding to 



- < 



(14) 



In terms of 7, this gives an upper bound of the value 
of the OAM in the plane of the crystal: 



< 27' 

and hence an image plane bandwidth 
i^ip = 1 + 472 . 



(15) 



(16) 



4.3 Far field bandwidth 

It is clear that in the far field of the plane of the crystal, 
instead of a real space argument, we need to use the angu- 
lar relationship P ^ a, expressed in IjlOp . where we apply 
the restriction for the maximum width of the detection 
modes given in (1141) : 



I I nkp 
'2L 



(17) 



From which, replacing Wg^i with Wp/^, we obtain an upper 
bound of the value of the OAM in the far field of the plane 
of the crystal: 



< 



27'ifl 



and therefore a far field bandwidth 

TT 



K 



FF 



1 



(18) 



(19) 



4.4 Measurement bandwidth 

If Kip and -f^FF are very different from each other, the 
resulting measurement bandwidth is given by the smaller 
of the two. For cases where the bandwidths are similar 
it is sensible to combine them. The convolution of two 
normal distributions of widths k and k' gives a normal 
distribution of width (fc~^ -I- fc'~^)~^/^. Similarly, we can 
get an estimate of the total measurement bandwidth by 
considering the convolution of two normal distributions 
of widths Kip and i^FF- The bandwidth of the resulting 
distribution is 



K 



1 + 47^ 



-1/2 



-2\ -1/2 



(20) 
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5 Analysis of the results 

For a comparison between the analytic and geometric ar- 
guments, we calculate the width of the distribution given 
by the modulus squared of the coefhcients in ^ and com- 
pare it to . In figure |3] we plot the two bandwidths 
as functions of Lr for 7 = 3, 7 = 5 and 7 = 7. The 
solid curves (red online) represent the measurement band- 
width calculated from the numerical evaluation of he an- 
alytical model. The dashed curves (green online) are the 
same bandwidths calculated with our geometrical argu- 
ment. The uppermost solid line (blue online) is the gen- 
eration bandwidth. Note that to achieve high dimensional 
entanglement the crystal length should be a small fraction 
of the Rayleigh range. 




0.000 0.002 0.004 0.006 0.008 0.010 



L 



R 



It is worth noting that for such value of 7 we have that 

H. There- 



Kip — Kpp 



— Kgcn, where i^gon is defined in 



fore in the optimal case we have K = Kg^n/V^- 

We define short crystal lengths as Lfj <^ 7r/47^, for 
which the generation bandwidth is large, meaning that 
the measurement bandwidth is dominated by the image 
plane overlap of the detection modes with the pump. This 
gives a measurement bandwidth of 



K « iCip = 1 + 472 



(22) 



Note that this short crystal limit is characterised by an 
independence of K on the crystal length. In fact, it can be 
seen in figure [3] that the leftmost part of the measurement 
bandwidth curves is flat (for the 7 = 7 curve this is not 
visible in this plot, but the slope of eq. ^IU\\ near the origin 
is zero for any 7), and that the range of values of Ln 
over which they stay flat is inversely proportional to 7^. 
For much longer crystals, ^ 7r/47'*, the measurement 
bandwidth, as modifled by the limiting overlap in the far 
field, becomes dominant, giving 



K « Kff = 1 



(23) 



In figure 21 we plot three different curves, that describe 
the value of the measurement bandwidth as a function 
of 7, for three different values of L/?. Note that for each 
choice of there is always an optimal value of 7 which 
maximises if, and it corresponds to the optimal value 
given in (PT|) . 



Fig. 3. (color online) The blue line (uppermost) is the 
generation bandwidth defined in (jl2p . the green curves 
(dashed) are calculated from our analytical treatment , 
and the red curves (solid) are the result of our geometrical 
argument. 

We see that the geometrical argument is in excellent 
agreement with the numerical evaluation of our analytical 
result. The effect of increasing 7 yields a higher measure- 
ment bandwidth for very small values of i^, but for large 
enough values of 7 and for fixed Lr, the measurement 
bandwidth eventually drops. Therefore it reaches a maxi- 
mum value for a particular crystal length. Under all con- 
ditions the measurement bandwidth never reaches that 
of the generation bandwidth, because we are restricting 
the measurement to modes with p — Q. Note, however, 
that the full generation bandwidth does not arise explic- 
itly from additional values of the 0AM but rather from 
entanglement in the radial quantum number p. 

Differentiation of eq. ([^0)1 with respect to the crystal 
length gives an estimate of the value of 7 corresponding 
to the highest measurement bandwidth for a given L^.- In 
this way we find 

7opt « (21) 



K 




2 4 6 8 10 12 ' 

Fig. 4. An example of a measurement bandwidth as a 
function of 7 for three different values of Lr: 0.001, 0.002 
and 0.003 (highest to lowest). 



It is not an easy matter to determine the requisite 
parameters for existing experiments. Most our own exper- 
iments, however, correspond to values of 7 in the range 
1.5 up to about 4. In order to achieve higher degrees of 
entanglement in 0AM, corresponding to larger Schmidt 
number, our analyses suggest that it would be desirable 
to press towards higher values of 7. 
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6 Conclusions 

We have shown two parameters determine the OAM band- 
width for entangled states produced by parametric down- 
conversion. Tliese parameters are the ratio of the widths 
of pump and detection modes 7 = Wp/wg^i, and the crystal 
thickness normalised to the Rayleigh range of the pump 

Lr = L/Zr- 

A simple geometrical argument approximates the ana- 
lytical results extremely well and allows us to suggest what 
needs to be adjusted in order to enhance the dimension- 
ality of the entanglement. We have restricted our analysis 
to a detection system that is sensitive to the LG p ~ 
modes only. It is for this reason that the measurement 
bandwidth can never reach that of the generation band- 
width for any combination of parameters. It is possible, 
however, to identify an optimum value of 7 to maximise 
the measurement bandwidth for any normalised crystal 
length Lr. 
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